Measurement of the Double Beta Decay Half-life of ^^"Te with the NEMO-3 Detector 



R. Arnold,^ C. Augier,^ J. Baker, ^'^^ A.S. Barabash,"* A. Basharina-Freshville,^ S.Blondel,^ M. Bongrand,^ 
G. Broudiii-Bay,6'7 y. Brudanin,^ A.J. Caffrey,^ A. Chapon,^ E. Chauveau/" D. Durand,^ V. Egorov,^ R. Flack,^ 
X. Garrido,2 J. Grozier,^ B. Guillon,^ Ph. Hubert,^'^ c.M. Jackson/" S. Jullian,^ M. Kauer,^ A. Klimenko,® 
O. Kochetov,® S.I. Konovalov,'* V. Kovalenko,^' D. Lalanne,^ T. Lamhamdi/^ K. Lang/^ Z. Liptak/^ 
G. Lutter,^''^ F. Mamedov," Ch. Marquet,^''^ J. Martin- Albo," F. Mauger,^ J. Mott,^ A. Nachabj^-"^ I. Nemchenok,^ 
C.H. Nguyen,^' 1^ F. Nova/^ P. Novella," H. Olisumi,!^ R.B. Palilka,!^ F. Perrot,^''^ F. Piquemal,^''^ J.L. Reyss,i^ 
B. Richards,^ J.S. Ricol,^''^ R. Saakyan,^ X. Sarazin,^ L. Siniard,^ F. Simkovic,!^ Yu. Shitov,^^!^ Smolnikov,^ 
S. S61dner-Rembold,i° I. Stekl,!^ J. Suhonen,^" C.S. Sutton,^! G. Szklarz,^ J. Tliomas,^ V. Tinikin,^ S. Torre,^ 
V.I. Tretyak,!'^ V. Uniatov,^ L. Vala," I. Vanyushin," V. Vasiliev,^ V. Vorobel,^^ Ts. Vylov,^'23 and A. Zukauskas^^ 

(The NEMO-3 Collaboration) 

^IPHC-DRS, Universite Louis Pasteur, CNRS, F-67037 Strasbourg, France 
^LAL, Universite Paris-Sud 11, CNRS/IN2P3, Orsay, France 
^INL, Idaho National Laboratory, Idaho Falls, Idaho 834-15, USA 
"^ITEP, Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 
^ University College London, London WCIE 6BT, United Kingdom 
^Universite de Bordeaux, CENBG, UMR 5797, F-33175 Gradignan, France 
CNRS/IN2P3, CENBG, UMR 5797, F-33175 Gradignan, France 
^ JINR, Joint Institute for Nuclear Research, 141980 Dubna, Russia 
^LPC, ENSICAEN, Universite de Caen, Caen, France 
^° University of Manchester, Manchester M13 9PL, UK 
" USMBA, Universite Sidi Mohamed Ben Abdellah, 30000 Fes, Morocco 
University of Texas at Austin, Austin, Texas 78712-0264, USA 
^^lEAP, Czech Technical University in Prague, CZ-12800 Prague, Czech Republic 
^'^IFIC, CSIC - Universidad de Valencia, Valencia, Spam 
Hanoi University of Science, Hanoi, Vietnam 
Saga University, Saga 840-8502, Japan 
^^LSCE, CNRS, F-91190 Gif-sur-Yvette, France 
^^Imperial College London, London SW7 2AZ, UK 
^^FMFI, Comenius University, SK-842 48 Bratislava, Slovakia 
Jyvdskyld University, 40351 Jyvdskyld, Finland 
'^^MHC, Mount Holyoke College, South Hadley, Massachusetts 01075, USA 
"^^ Charles University in Prague, Faculty of Mathematics and Physics, CZ-12116 Prague, Czech Republic 

^'^ Deceased 
(Dated: January 19, 2013) 

We report results from the NEMO-3 experiment based on an exposure of 1275 days with 661 g 
of ^^"Te in the form of enriched and natural tellurium foils. The double beta decay rate of ^"^"Te 
is found to be greater than zero with a significance of 7.7 standard deviations and the half-hfe is 
measured to be Ty2 — (7-0 ± 0.9(stat) ± l.l(syst)) x 10^" yr. This represents the most precise 
measurement of this half-life yet published and the first real-time observation of this decay. 

PACS numbers: 23.40.-s, 21.10.Tg, 14.60.Pq 



The first evidence of double beta decay (2i//3/3) ap- 
peared in 1950 through the observation of ^"^"Xe from 
the decay of ^'^"Te in rock samples [1]. This result 
was met with scepticism for the ensuing 15 years un- 
til the results of a number of other geochemical exper- 
iments began to confirm the observation. There was, 
however, significant disagreement between two distinct 
sets of these measurements that was not immediately 
resolved. Several groups measured a long half-life of 
ss 2.7 X 10^^ yr [2, 3] while others obtained a significantly 
shorter half-hfe of « 0.8 x 10^^ yr [4-7]. One hypothe- 
sis to explain the difference is based on the observation 
that shorter half-lives were measured in rock of relatively 



young age (~ 10^ — 10^ yr), while the longer half-lives 
were measured in r elatively old rock (~ 10^ yr) [2, 4]. It 
has even been suggested that there is a time dependence 
in the value of the weak interaction coupling constant [8] . 
Recent papers [9] attempt to explain the long-held dis- 
crepancy between these measurements as being caused by 
catastrophic xenon loss in the older samples. To date, the 
only direct evidence for the ^■^"Te 2z^/3/3 process in an ex- 
periment comes from MIBETA, which reported a half-life 
of (6.1 ± 1.4(stat)t3j(syst)) x 10^° yr [10] by comparing 
different crystals isotopically enriched in ^-^"Te and ^^^Te, 
assuming that any difference in rate was due to 2vPj3 
events (^^^Te has a much longer 2zy/3/3 half-life). How- 
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ever, a systematic uncertainty of about 50% rendered this 
result somewhat inconclusive. In this Letter, we present 
the first direct, high-precision measurement of 2vPP de- 
cay of i30Te with the NEMO-3 detector. In addition, a 
search for neutrinoless double beta decay {Oi^/3/3) and for 
the decay with Majoron emission {Qv^j3j3) is reported. 

The NEMO-3 detector is located in the Modane Un- 
derground Laboratory. The detector [11] contains almost 
9 kg of seven different /3/3 isotopes in the form of thin foils. 
It provides direct detection of electrons from the double 
beta decay by the use of a tracking device based on open 
Geiger drift cells and a calorimeter made of plastic scin- 
tillator blocks coupled to low-radioactive photomultipli- 
ers (PMTs). For 1 MeV electrons the timing resolution 
is 250 ps and the energy resolution (full width at half 
maximum) is about 15%. A magnetic field surround- 
ing the detector provides identification of electrons by 
the curvature of their tracks. In addition to the electron 
and photon identification through tracking and calorime- 
try, the calorimeter measures the energy and the arrival 
time of these particles while the tracking chamber can 
measure the time of delayed tracks associated with the 
initial event for up to 700 /iS. The calorimeter energy 
scale is calibrated approximately once per month using a 
^'^^Bi source providing conversion electrons of 482 keV 
and 976 keV (K-lines). Stability of the calorimeter re- 
sponse is surveyed twice a day by a laser system. The 
data presented in this Letter correspond to 1275 days of 
data taking between October 2004 and December 2009. 
Two different foils are used in the analysis: a Te foil, en- 
riched at a level of 89.4 ± 0.5% corresponding to 454 g of 
"OTe, and a natural Te foil which contains 33.8% ^^OTe, 
corresponding to 207 g of ^■^'^Te. 

When searching for rare processes, the background es- 
timation is paramount as it will limit the final sensitiv- 
ity of the experiment. An exhaustive program of work 
has been carried out to measure the very large number 
of sources of background present in the NEMO-3 detec- 
tor. The method of the background measurement and 
its validation with a highly radiopure Cu foil is described 
in [12]. There are three categories of background: the 
external background, originating from radioactivity out- 
side the tracking chamber; the tracking volume back- 
ground, which includes radon in the tracking gas and 
the drift cell wire contamination; and the internal back- 
ground due to radioactive impurities inside the source 
foils whose dominant isotopes are ^°K, ^'^'*™Pa, ^^°Bi, 
^^''^Bi and ^"^Tl. The different background contributions 
are estimated by measuring independent event topolo- 
gies, both for enriched and natural Te. 

The external background originates from components 
outside the tracking volume. PMTs are the main con- 
tributors since they have glass and electronic components 
that, though at a very low level, contain traces of ^^'^Bi, 
^"^^Tl and ''^K. The external background is measured 
with 7-ray Compton scattering in the scintillators, either 



producing an electron that crosses the tracking chamber 
or depositing energy in one scintillator, followed by an 
electron emitted from the source to another scintillator. 
The reliability of the external background model is il- 
lustrated by the energy distributions of the one-electron 
crossing events in Figs. 1(a) and (b). 




FIG. 1. (Color online) Energy sum distribution for 
crossme- electron events for (a) ^^°Te and (b) "="Te (b); 
Energy distribution of electron events coming from the 
source foil for (c) "°Te and (d) "''^Te. Dots correspond 
to the data and histograms to the fit of the background 
model. 



The source foil activities in ^'^^'"Pa, '^''K and ^^°Bi are 
determined with single-electron events coming from the 
foil. The energy distribution of the observed events and 
the result of the fit of the different components of the 
background are presented in Figs. 1(c) and (d). The good 
agreement between the data and the fit demonstrates the 
reliability of the internal background model. The foil 
activity in ^^^Bi is measured using events with a single 
electron accompanied by a delayed a-particle track. This 
topology is a signature of the /3 decay of ^^^Bi to ^^"^Po 
followed by a decay of ^^'^Po to ^^°Pb. The foil activity 
in 208 rp]^ measured with events that contain one elec- 
tron and either two or three photons emitted from the 
foil. The results of measurements of the internal con- 
tamination by ^^'^Bi and ^os^j g^j.g reported in [12]. The 
measured foil activities are summarized in Table I. 

The background from single /3 decay in the tracking 
volume is of importance if the decay occurs near the foil. 
The main source of this background is due to daughters 
of radon: ^^'^Pb, ^^"^Bi and ^^°Bi. The radon activity in 
the tracking chamber is measured using ea events as for 
the measurement of internal ^^^Bi background. The dis- 
tribution of these events is measured as a function of the 
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Impurity 








0.17 ± 0.04 


0.24 ± 0.04 




0.29 ±0.05 


0.25 ±0.12 


234m 


2.49 ± 0.05 


0.73 ±0.04 




19.9 ±0.4 


18.4 ±0.3 


«K 


14.7 ±0.2 


12.3 ±0.2 



TABLE I. Background contaminations measured in the Te 
source foils (in mBq). 



location in the tracking volume. For the data presented 
here, the mean ^^^Rn activity in the whole gas volume 
is 209 ± 2 mBq. Unlike the preceding radon daughters, 
^^"Pb has a long half-life of about 22 years. It is therefore 
not in equilibrium with ^^^Rn in the tracking volume and 
most ^^"Pb was deposited during the construction of the 
detector. ^^°Bi produced in /3 decay of ^^"Pb contributes 
to the low-energy background below 1 MeV. The ^^"Pb 
deposition on drift cell wires, measured by detecting elec- 
trons from ^^"^Bi f3 decay, was found to be vary signifi- 
cantly in different sectors [12]. In contrast to ^^^Rn the 
concentration of ^^°Rn in NEMO-3 is very small and its 
contribution to the total background in the ^■^^Te sectors 
is less than 1%. 

The measured activities are used to estimate the back- 
ground contribution in the two-electron channel with a 
Monte Carlo (MC) simulation. Signal and background 
MC events are generated using a GEANT-based simula- 
tion [13] of the detector with the initial kinematics given 
by the event generator decayO [14]. 

The two-electron events are selected with the follow- 
ing requirements. Two tracks of a length greater that 
50 cm with curvature corresponding to a negative charge 
are reconstructed. Both tracks originate from a common 
vertex in the foil and terminate in isolated scintillators 
with a single energy deposit greater than 0.2 MeV. The 
time-of-flight information is consistent with the hypoth- 
esis that two electrons were emitted from the same point 
on the source foil. No photon or delayed a track is de- 
tected in the event. These selection criteria lead to a 
2^/3(3 detection efficiency of 3.5% for enriched ^'^'^Te and 
2.8% for "'^'Te, where the difference is mainly due to the 
source foil thickness. Several factors contribute to this 
low efficiency, but the most important are the geomet- 
rical acceptance of the detector, the effect of the energy 
threshold and the tracking algorithm inefficiencies. 

We measure the ^'^°Te half-life with data from the 
enriched Te source foil by performing a likelihood fit 
to the binned energy sum distribution in the interval 
[0.9 — 2] MeV. This interval is chosen using the MC simu- 
lation to maximize the signal significance. It reduces the 
2i'/3f3 efficiency by a factor of 0.7. The result of the fit is 
presented in Figs. 2(a), (c), and (e). Using the MC sim- 
ulation the number of background events in the interval 
[0.9 — 2] MeV is estimated to be 363 ± 25 events, of which 
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FIG. 2. (Color online) (a,b) Distribution of the sum 
of the electron energies; (c,d) individual electron en- 
ergy; and (e,f) cosine of the angle between the two elec- 
tron tracks for two-electron events selected from the two 
Te foils: (a,c,e) enriched in ^^°Te and (b,d,f) natural Te. 



141 events are associated with the external background, 
179 with the internal background and 43 with radon in- 
duced background. The number of events in excess of the 
background in the interval [0.9 — 2] MeV is determined 
to be 



n{2iyl3l3) = 178 ± 23, 



(1) 



with a signal significance of 7.7 standard deviations and 
a signal-to-background ratio of S/B = 0.5. This corre- 
sponds to a "OTe half-hfe of 



T^J.^ ^ (7.0 ± 0.9) X 10^" 



(2) 



The main systematic uncertainty on the measured 
^^°Te half-life is associated with the background esti- 
mation and is due to the small signal-to-background ra- 
tio. The uncertainty on the number of expected back- 
ground events has been obtained by applying the largest 
variations of the component activities in the background 
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model. The corresponding uncertainty on the -'^^''Te half- 
life is 14%. Another systematic uncertainty is associ- 
ated with the two-electron detection efficiency in NEMO- 
3 which is found to be correct within an accuracy of 5%. 
This uncertainty is determined with a calibrated ^"''Bi 
source and a dedicated ^°Sr source which decays to ^"^Y, 
a pure fi emitter of Qp = 2.28 McV. Finally, the source 
foil thickness and the GEANT model of electron energy 
losses in dense thin media contribute a systematic un- 
certainty of 4% which is estimated by comparing signals 
from metallic and composite ^"^Mo source foils [15]. The 
total systematic uncertainty of 15% is obtained by adding 
the individual contributions in quadrature. 

The measurement of the /3/32i/ half-life is verified using 
the natural Te foil. The energy and angular distribu- 
tions of the two-electron events of the natural Tc data 
are presented in Fig. 2b,d,f and are compared with the 
expected MC distribution using the measured half-life 
T^J^ = 7.0 X 10^° yr and the background model of the 
natural Te foil. There are 65 ± 8 /3/3 events and 316 ± 28 
background events expected in the electron energy sum 
interval [0.9 - 2] MeV. The total number of 381 ± 29 ex- 
pected events is in good agreement with the 377 observed 
events. 

The '^'^"Te data (Fig 2a) are also used to set a 
limit on the Oz^/3/3 and Ovx^PP processes with the CL^ 
method [16]. The method uses the full information of 
the binned energy sum distribution for signal and back- 
ground, as well as the statistical and systematic uncer- 
tainties and their correlations as described in [17]. 

The total efficiency to detect Oj^/3/3 decay of ^^°Te is 
estimated to be (13.9 ± 0.7)% yielding a limit of 

T°/2 > 1.3 X 10^3 yr (90% C.L.), (3) 

which is an order of magnitude less stringent than the 
Umit obtained by the CUORICINO Experiment [18] 
based on 11 kg of ""Te. 

The detection efficiency for the decay with ordinary 
(spectral index n = 1) Majoron emission (see discussion 
in [19] and references therein) is (9.6 ± 0.5)% and the 
limit is determined to be 

T^J^" > 1.6 X 10^2 yr (90% C.L.), (4) 

which is a factor of 7 more stringent than the previous 
best limit from MIBETA [10]. The corresponding limit 
on the coupling constant of the Majoron to the neutrino 
is goo < (0.6 — 1.6) X 10^^ (using nuclear matrix elements 
from [20-24]) and is comparable with the best present 
limits. 

In summary, the 2z//3/3 decay ^'^"Te half-life measured 
with the NEMO-3 detector is 

T^'f^ = (7.0 ± 0.9 (stat) ± 1.1 (syst)) x 10^° yr. (5) 



With this result, the corresponding nuclear matrix ele- 
ment can be extracted according to 

(1^1%)"' =G^"|M2"^ (6) 
where C^" = 4.8 x 10" yr^i (for qa = 1.254) is the 
known phase space factor [25], which yields the result 
M"^" = 0.017 ±0.002 (scaled by the electron mass). This 
value for M^'' may be used to fix the gpp parameter of 
the QRPA theory, which corresponds to the strength of 
the nucleon-nucleon interaction inside the nucleus. It has 
been suggested that this will improve the M*^" calcula- 
tions [20, 21, 26]. 

The NEMO-3 result for the ^^°Te half-life is consis- 
tent with the geological measurements made in younger 
rock samples and is the most precise measurement of this 
isotope half-life to date. 
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